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Abstract 



The classical problem of the thermal explosion in a long cylindrical vessel 
is modified so that only a fraction a of its wall is ideally thermally conduct- 
ing while the remaining fraction 1 — a is thermally isolated. Partial isolation 
of the wall naturally reduces the critical radius of the vessel. Most inter- 
esting is the case when the structure of the boundary is a periodic one, so 
that the alternating conductive a and isolated 1 — a parts of the boundary 
occupy together the segments 2ir/N (N is the number of segments) of the 
boundary. A numerical investigation is performed. It is shown that at small 
a and large N the critical radius obeys a scaling law with the coefficients 
depending upon N. For large N is obtained that in the central core of the 
vessel the temperature distribution is axisymmetric. In the boundary layer 
near the wall having the thickness m 2-irro/N (ro —the radius of the vessel) 
the temperature distribution varies sharply in the peripheral direction. The 
temperature distribution in the axisymmetric core at the critical value of 
the vessel radius is subcritical. 
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1 Introduction 



We revisit in this Note the classical problem of a thermal explosion in a long circular cylin- 
drical vessel containing an exothermically reacting gas at rest. It was formulated and solved 
under some assumptions by D. A. Frank-Kamenetsky Jl|, (see also ||). In the original 
formulation it was assumed that the wall of the vessel is ideally conducting, so that the gas 
temperature at the boundary is equal to the temperature of the ambient medium. Such 
boundary condition made the problem cylindrically symmetric, and the symmetry essen- 
tially simplified the solution. In the present Note the problem is modified in the following 
way: the wall is partially isolated so that the symmetry is lost. The critical values of the 
radius of the vessel are determined numerically. Especially instructive results were obtained 
for the cases when the isolated parts of the boundary are distributed periodically with large 
angular frequency, and the isolated part of the boundary is large. The scaling laws for the 
critical values were found. For large angular frequencies it was found that there exists an 
axisymmetric core of the temperature distribution which occupies a major part of the vessel. 
The conditions in this core at the critical case were found to be subcritical. 



2 Mathematical Problem Formulation 

Assume that a gas at rest is enclosed in a long cylindrical vessel of radius tq. An exothermic 
reaction is going in the gas with the thermal effect Q per unit mass of reacted gas. For the 
reaction rate the Arrhenius law is assumed with the activation energy E. If the thermal 
effect and the activation energy are large, it can be shown (see ||) that an 'intermediate- 
asymptotic' steady state regime is achieved. For this regime the gas consumption in the 
reaction can be neglected, and, from the other side, the temperature distribution in the 
vessel is steady. Applying the Frank-Kamenetsky large activation energy approximation, a 
non-linear equation for dimensionless reduced temperature u is obtained: 

Am + AV = 0, (2.1) 

where 

(T - To)E 

the Laplace operator A is related to dimensionless variables p = r/r , 8; r,8 are the polar 
coordinates. The constant A is 

A = j, 1= (e^KRT%c/QEa(T Q )} 112 . (2.3) 

Here T is the absolute temperature, To-the temperature of the ambient medium, i?-the 
universal gas constant, K-the molecular thermal diffusivity, c-heat capacity of gas per unit 
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volume, o~(T) is the pre-exponential factor in the Arrhenius reaction rate expression: a slow 
function of temperature. 

In the classical problem formulation it was assumed that the whole wall of the vessel 
is ideally heat conducting, so that the gas temperature at the boundary is equal to the 
temperature of the ambient medium. This gives a Dirichlet condition for the equation ( |2.1| ): 

u(l,0) = O. (2.4) 

The boundary value problem formulation under this condition is axisymmetric, and 
this was important for obtaining the analytic solution in an explicit form. D. A. Frank- 
Kamenetsky showed that the solution to the problem fl2.1|) , (|2.4j) does exist for A < A cr = \[2 
only. Physically it means that a quiet, non-explosive proceeding of the reaction is possible 
only if the radius of the vessel is less than a critical one: tq < (ro) cr = \[2l. This condition 
is known as the condition of the thermal explosion. 

In the present Note the following modification of the problem ( |2.1|) , (2.4) is proposed: 



Only a fraction a of the wall is heat conducting, while the fraction 1 — a is thermally isolated. 
The simplest formulation corresponds to a mixed problem (Figure l,a) 

u(l,0) = O at < 9 < 27ra . . 

d p u{l,6) = at 2na<9<2n. v ' ; 

More interesting is the case when the isolated part of the wall is not concentrated on a 
single arc, but is distributed periodically (Figure l,b): the boundary p = 1 is divided into 
N segments 

2tt 2tt 

^(l + k)>9>—k, fc = 0,l,...,iV-l. (2.6) 

The fraction a of each segment is left heat conducting, while the fraction 1 — a becomes 
isolated. In this case the mixed boundary condition at p = 1 has the form: 



„s 27T , . 2ira 2n , 

u (i,e) = a at —k<e< — + —k, 

d p u(l,0) = O at 2 -^ + ^k<6<^(k + l) 



(2.7) 



The central question addressed in the present Note is: What are the asymptotic laws for 
the critical radius if A^ — > oo and a — > 0, i.e. the period of the boundary condition tends to 
zero, and the isolation is close to a complete one. Remember: for the complete isolation the 
critical radius is equal to zero. 
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3 The Numerical Method 



In order to answer the questions posed above, we must numerically evaluate the critical 
values X cr for fixed N and a. There are two special items in determining A cr : the singularity 
of the linearized equation (|2.1| ) : 

A5u + X 2 e u 5u = (3.1) 

at the critical value of A = A cr , and the dependence of the values of A cr obtained by dis- 
cretization upon the number of points, n used to discretize each dimension. Here 5u is the 
perturbation of the solution. 

For fixed values of N, n, and a, we determined a trajectory of solutions u versus A by 



solving equation (2.1) with a Newton-Raphson method. This method requires the solution 



of the linearized equation (|3.1| ) which becomes singular at A = A cr . In practice it prevents 
us from approaching the critical point closely. Therefore to make an accurate determination 
of A cr an extrapolation procedure was used. It was assumed that for A approaching A cr a 
parabolic approximation is valid: 

A c 2 r -A 2 = C(H- Mo ) 2 . (3.2) 

The parameters A cr , C and Uq were determined to fit the last 10 points on the trajectory u 
versus A approaching A cr . The approximation ( |3.2|) happened to be satisfactory. Typically 
the fit ( p.2| ) is accurate to a few parts in 10 6 . 

The above procedure yielded an estimate for A cr as a function of N, a and the number of 
discretization points n. It is natural to remove the dependence of A cr of non-physical, com- 
putational parameter n. For this purpose another extrapolation was used. In our numerical 
approximations the second-order accurate discretizations of the operators was employed. If 
the boundary conditions were smooth, the solution u would approach a limit with an error 
of the order of 1/n 2 . But the boundary conditions are non-smooth, the derivative of the 
solution is discontinuous at the boundary, and it causes the order of the approximation to 
decrease. Extensive numerical calculations have shown \ 2 cr to vary linearly with 1/n. There- 
fore, the following iterative procedure was used: the value A cr (iV, a, n) in ( |3.2j ) was calculated 
for three different values of n. These three values are then fit to a linear function a + b/n. If 
the fit was poor, as might happen if the values of n were too small, the procedure is repeated 
with larger values of n and so on, until a satisfactory fitting was obtained. The extrapolation 
n —>■ oo is simply the value of the coefficient a. 



4 Results of the Numerical Analysis 

The results obtained by numerical solution are represented on Figures 2-4. Three instructive 
properties are revealed. 
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(i) On the graph of Figure 2 the values of \ 2 cr are presented for growing values of A" as the 
functions of 1/a. It is seen that the critical value X cr , i.e. the critical radius of the vessel for 
large N is practically insensitive to a up to a very small. For small A" the dependence of A cr 
on a is strong. Clearly, for any N, X cr = for a = 0, but it is instructive that for instance, 
for N = 256 when only 1/512 (0.2 percent) part of the boundary is heat conducting, the 
critical value of the radius is only 4 percent less than the critical radius for wholly heat 
conducting wall. 

(ii) For large N, starting, say, from iV = 32, there exists an internal core < p < p* 
where the solution is close to axisymmetric one (see Figure 3). The value p* was selected so 
that |w m ax(p*? 0) — Wmin(p*, 6) \ < 10~ 4 . Introducing the mean value = (ti mffl + w m i n )/2 we 
notice, that for < p < p* the solution is close to axisymmetric, so that the equation (|2.1| ) 
and the boundary condition at p = p* assume the form: 

1 d du o 

z-p — + A e = 0, u = at p = p*. (4.1) 

p dp dp 



Transforming the variables 



u = Ujf + Vy R=P- (4.2) 



p* 



we reduce the problem to a classic one: 

where < R < 1, A 2 = X 2 ple u *. We calculate now the values of A 2 r = \\ r p\t u * where A cr is 
the critical value obtained in previous calculations. The graphs A 2 r as the functions of 1/a 
for different are presented on Figure 4. It can be seen that the values of K 2 cr are always less 
than 2 (within the limits of our numerical accuracy). This means that the 'rugged' boundary 
layer near the wall p = 1 controls the approach to criticality. The thickness of this layer is 
of the order of the length of the segment 2ir/N . The angular derivative dgu in the boundary 
layer is large. The value decreases as A^ increases. 

(iii) The intermediate power laws are observed for \ 2 cr at large A^ and small a: 

\ 2 cr = S{N)a tiN) . (4.4) 

The values of S(N) and t(N) for various A^ are given in the Table. 

5 Conclusion 

Non-axisymmetric modification of the problem of thermal explosion in a cylindrical vessel 
is formulated. The boundary is partly isolated, and only partly ideally conducting. Special 
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attention is paid to the case of periodic distribution of the isolated and conducting parts. 
The critical values of radius and other relevant properties are obtained numerically. 

It is shown that for the period small in comparison with the vessel radius the critical 
value of radius of the vessel is practically insensitive to the relative size of the open area 
of the wall up to its very small values. The temperature distribution in the central core 
is axisymmetric and subcritical even at globally critical conditions: the criticality is due 
to a thin boundary layer near the wall where the temperature distribution is highly non- 
axisymmetric. Intermediate power laws are obtained for the critical radius as the function 
of the relative open area of the wall. 

This work was supported in part by the Applied Mathematical Sciences subprogram of 
the Office of Energy Research, U.S. Department of Energy under Contract DE-AC-03-76- 
SF00098, and in part by the National Science Foundation under Grants DMS94-14631 and 
DMS-2732710. 
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Table 

N S(N) t(N) 

32 2.03 0.10 

64 2.03 0.055 

128 2.02 0.028 

256 2.01 0.015 
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Figure Captions 



Figure 1. A fraction of the wall is isolated, (a) Isolated (5/6) and conducting (1/6) parts 
are connected, (b) Isolated and conducting parts are distributed periodically. 

Figure 2. Dimensionless critical radius as the function of conducting fraction a for different 
angular frequencies. It is seen that at large frequencies the critical radius is practically 
a-independent up to very small values of the conducting fraction a. 

Figure 3. The solution reveals an axisymmetric internal central core and 'rugged' boundary 
layer (TV = 32, a = 1/32). 

Figure 4. The temperature distribution in the internal central core is a subcritical one: 
A 2 < 2. 
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